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THE DOMAIN OF ANALYTICITY OF SOLUTIONS TO THE THREE-DIMENSIONAL 

EULER EQUATIONS IN A HALF SPACE 

IGOR KUKAVICA AND VLAD VICOL 



O 

Abstract. We address the problem of analyticity up to the boundary of solutions to the Euler equations 
in the half space. We characterize the rate of decay of the real-analyticity radius of the solution u(t) in 
terms of exp L || Vu(s)||£oo ds, improving the previously known results. We also prove the persistence of the 
sub-analytic Gevrey-class regularity for the Euler equations in a half space, and obtain an explicit rate of 

£f) ' decay of the radius of Gevrey-class regularity. 

1. Introduction 

p~H The Euler equations on a half space for the velocity vector field u(x, t) and the scalar pressure field p(x, t), 

where x G Vt = {x E R 3 : X3 > 0} and t > 0, are given by 

d t u + (u ■ V)u + Vp = 0, in Q x (0, 00), (E.l) 

V-u = 0, in ft x (0,oo), (E.2) 

u-n = 0, on dVL x (0,oo), (E.3) 

where n = (0,0,— 1) is the outward unit normal to dfl = {x <E M 3 : X3 — 0}. We consider the initial value 
problem associated to (|E.1[) - (|E.3[) with a divergence free initial datum 

u(0) = mo, in il. (E.4) 

The local existence and uniqueness of 7? r -solutions, with r > 3/2 + 1, on a maximal time interval [0,T„) 
holds (cf. OH [HI [23J E2 [36]), and lim T ^T. Jq II cwlu(t)\\ L °°dt = 00, if T* < 00 (cf. [6]); additionally the 
f — ■ persistence of C°° smoothness was proven by Foias, Frisch, and Temam [17]. In this paper we address the 

solutions of the Euler initial value problem evolving from real-analytic and Gevrey-class initial datum (up to 
the boundary), and characterize the domain of analyticity. We emphasize that the radius of real-analyticity 
gives an estimate on the minimal scale in the flow [221 [3S] , and it also gives the explicit rate of exponential 
decay of its Fourier coefficients [18] . 

In a three dimensional bounded domain, the persistence of analyticity was proven by Bardos and Bc- 
nachour [3] by an implicit argument (see also Alinhac and Metivier [I]). In [H [7] the authors give an 
explicit estimate on the radius of analyticity, but which vanishes in finite time (independent of T*). How- 
ever, the proof of persistency [2] can be modified to show that the radius of analyticity decays at a rate 
proportional to the exponential of a high Sobolev norm of the solution (see also \T). On the three dimen- 
sional periodic domain (or equivalently on R 3 ) this is the same rate obtained by Levermore and Oliver in 
[30j , using the method of Gevrey-class regularity. This Fourier based method was introduced by Foias and 
Temam [18] to study the analyticity of the Navier- Stokes equations. For further results on analyticity see 
[3[Il[Tl[ini[lTl[21[251[151[221[2Zl[S3[SS]. Explicit and even algebraic lower bounds for the radius of 
analyticity for dispersive equations were obtained by Bona, Grujic, and Kalisch in [9l[10] (see also [8lfTl]V 

In [2f[ we have proven that in the periodic setting, or on E 3 , the analyticity radius decays algebraically 
in the Sobolev norm || curlu(£)||ijr-, with r > 7/2, and exponentially in J Q ||Vu(s)||l°c, for all t < T*. In 
the present paper we show that the algebraic dependence on the Sobolev norm holds in the case when the 
domain has boundaries (cf. Theorem 12. ip . thereby improving the previously known results. The interior 
analyticity in the case of the half-space, for short time (independent of T*), was treated in [35]. We note 
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that the shear flow example of Bardos and Titi [5] (cf. [TU [35]) m &y be used to construct explicit solutions 
to the three-dimensional Euler equations whose radius of analyticity is decaying for all time. 

Additionally we prove the persistence of sub-analytic Gevrey-class regularity up to the boundary (cf. p~8j 
[31] ) for the Euler equations on the half space. To the best of our knowledge this was only known for the 
periodic domain cf. [26j[30], but not for a domain with boundary. The methods of [IJ [3j [27l [35] rely essentially 
on the special structure of the complex holomorphic functions, and do not apply to the non-analytic Gevrey- 
class setting. 

The presence of the boundary creates several difficulties that do not arise in the periodic setting. In 
particular we cannot use Fourier-based methods, nor can we use the vorticity formulation of the equations. 
Instead we need to estimate the pressure, which satisfies (cf. [31)) the elliptic Neumann problem 

— Ap = djUidiUj, in fi x (0,oo), (P-l) 

-£ = (u-V)wn = 0, onSQx (0,oo), (P.2) 

an 



since n — (0,0,-1), where the summation convention on repeated indices was used in (|P.1|) . In order to 
close our argument we need to show that the pressure has the same analyticity radius as the velocity, and 
so we cannot appeal to the inductive argument of Lions and Magenes |31) . Moreover, the nature of the 
elliptic/hyperbolic boundary value problem imposes certain restrictions on the weights of the Sobolev norms 
that comprise the analytic norm. The analytic norm we define (cf. Section [3]) respects the symmetries of the 
problem and is adequate to account for the transfer of derivatives arising in the higher regularity estimates 
for the pressure. 

The paper is organized as follows. In Section [5] we state our main result, Theorem 12.11 In Section [3] 
we prove the main theorem assuming two key estimates on the convection term and the pressure term, 
Lemma l3.1l and Lemma l3.2l Section [4] contains the proof of the commutator estimate Lemma l3.ll and lastly, 
the higher regularity estimates for the pressure and the proof of Lemma l3.2l are given in Section [5] 

2. Main Theorems 

The following statement is our main theorem addressing the analyticity of the solution. Theorem 12.21 
below concerns the Gevrey class persistence. 

Theorem 2.1. Fix r > 9/2. Let uq € H r (Q) be divergence-free and uniformly real-analytic in Q. Then 
the unique solution u{t) £ C(0,T*; H r (fl)) of the initial value problem associated to the Euler equations 
(|E.1I) - (|E.4I) is real-analytic for all time t < T*, where T* G (0, oo] denotes the maximal time of existence of 
the H r -solution. Moreover, the uniform radius of space analyticity r(t) of u(t) satisfies 

r(t) > — L-- exp (-C J HVutolli-ds) , (1) 

where C > is a constant that depends only on r, while Cq has additional dependence on uq as described in 
([TTj) below. 

Remark 1. The lower bound ([1]) improves the rate of decay from Bardos and Benachour 3 on a bounded 
domain (which can be inferred to be proportional to exp/ Q ||w(s)||jjr<is), and it matches the rate of decay 
we obtained in [26] on the periodic domain. 

Remark 2. The proof of Theorem 12 . 1 1 also works in the case of the half-plane (recall that in two dimensions 
T* may be taken arbitrarily large, cf. |32[|37j ) with the same lower bound on the radius of analyticity of the 
solution. Since in two dimensions ||Vit(£)||i,°o grows at a rate of Cexp(Ct), for some positive constant C, the 
estimate ([IJ shows that the rate of decay of the analyticity radius is at least Cexp(— Cexp(Ci)), for some 
C > 0. This recovers the two-dimensional rate of decay obtained by Bardos, Benachour and Zerner [4 on a 
bounded domain and by the authors of this paper on the torus [26J . It would be interesting if one could prove 
a similar lower bound to ([IJ but where the quantity J ||Vu(s)||l°o ds is replaced by J Q || curlw(s)||L°o ds. In 
particular, such an estimate would imply in two dimensions that the radius of analyticity decays as a single 
exponential in time. 
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Recall (cf. [3T]) that a smooth function v is uniformly of Gevrey-class s, with s > 1, if there exist M, t > 
such that 

\d<*v(x)\ < M^, (2) 

for all x € 17 and all multi- indices aeNj. When s = lwe recover the class of real-analytic functions, and 
for s € (1, oo) these functions are C°° smooth but might not be analytic. We call the constant r in $2§ the 
radius of Gevrey-class regularity. The following theorem shows the persistence of the Gevrey-class regularity 
for the Euler equations in a half-space. 

Theorem 2.2. Fix r > 9/2. Let uq be uniformly of Gevrey-class s on 17, with s > 1, and divergence-free. 
Then the unique H r -solution u(t) of the initial value problem (1E.1JI - (1E.4|I on [0,T„) is of Gevrey-class s, for 
all t < T*, and the radius t(£) of Gevrey-class regularity of the solution satisfies the lower bound (flj). 

3. Proofs of Theorem 12. II and Theorem 12.21 

For a multi-index a — (0:1,02,03) in Nq, we denote a' — (01,02). Define the Sobolev and Lipshitz 
semi-norms | • | m and | • | mj00 by 

\v\ m = Y, M a \\d a v\\ L 2, (3) 

\a\—m 

and 

|«koo= J2 M a \\d a v\\ L ~, 

|a|— m 

where 

v'M /oi + «2 



a'\ \ ai / 

The need for the binomial weights M a in ([3]) shall be evident in Section[5]where we study the higher regularity 
estimates associated with the Neumann problem (|P.1[) - (|P.2[) . For s > 1 and r > 0, define the space 



X T = {u G C°°(fi) : ||v|| Xt < 00}, 
where 

00 r m-3 

HI** = E ^(^F 

m— 3 v ' 

Similarly let F T = {u e C°°(17) : |H|y T < 00}, where 

2, (m-3)r m - 4 



My t 



£*>> ' (m-3)! s 

m=4 v ' 



Remark 3. The above defined spaces X T and Y T can be identified with the classical Gevrey-s classes 
as defined in 3f . On the full space or on the torus, the Gevrey-s classes can also be identified with 
Z>((-A) r / a exp(r(-A) 1 / 2 »)) ( c f. |fg| |2g| gjg ) . 



We shall prove Theorems 12. fl and 12.21 simultaneously by looking at the evolution equation in Gevrey-s 
classes with s > 1. If uq is of Gevrey-class s in 17, with s > 1, then there exists r(0) > such that uo € -Xt(o), 
and moreover t(0) can be chosen arbitrarily close to the uniform real-analyticity radius of uo, respectively 
to the radius of Gevrey-class regularity. Let u(t) be the classical 7? r -solution of the initial value problem 

With the notations of Section [5] we have an a priori estimate 



l,\Ht)\\ XHt) -r(t)||«(t)|| n(t) + £ (|l«(t)U) ^^7. (5) 

m=3 V ' V '' 
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Fix m > 3. In order to estimate (d/dt)\u(t)\ m , for each |a| = m we apply d a on (jE.ip and take the L 2 -inner 
product with d a u. We obtain 

\j\\d a u\\ 2 L2 + < d a (u ■ \7u),d a u > + < Vd a p,d a u >= 0. (6) 

On the second term on the left, we apply the Leibniz rule and recall that < u ■ V9 Q w, d a u >— 0. For the 
third term on the left of (|6]) we note that since n = (0, 0, —1) and u ■ n = on dfl, we have that d a u ■ n = 
for all a such that a 3 — 0. Together with V ■ u — in $7 this implies that < Wd a p,d a u >— whenever 
a 3 = 0. Using the Cauchy-Schwarz inequality and summing over \a\ = m we then obtain 

j t Wm< ]T ]T M a ( a \\d^u-wd a -^u\\ L ,+ Yl McWWpU*. 

\a\=m p<a,Pj^Q ^ P ' \ a \=m,a a ^Q 

Combined with ([5]), the above estimate shows that 

^||«(*)||x tW < r(«)||«(t)||y TW + C + V, (7) 

where the upper bound on the commutator term is given by 

oo / \ m — 3 

m=3|a|=m^<a,/9/0 VK/ v ; 

and the upper bound on the pressure term is 

00 r m-3 

m=3 |a|=m,a 3 ^0 

In order to estimate C we use the following lemma, the proof of which is given in Section [4] below. 
Lemma 3.1. There exists a sufficiently large constant C > such that 

C<C(Ci+C 2 \\u\\ Yr ), 

where 

C\ = Ml,ooM 3 + M2,ooM2 +t|u| 2 ,ooM3, 

and 

C 2 = r|u|i >00 + t 2 \u\ 2 ,oo + t 3 M 3 ,oo + r 3/2 ||u||x r - 
The following lemma shall be used to estimate V. The proof is given in Section [5] below. 
Lemma 3.2. There exists a sufficiently large constant C > such that 

V<C(V 1 +V 2 \\u\\ Yt ), 
where 

V\ = Ml,ooM3 + |«|2,ooM2 + t|w| 2 ,oo|"|3 + 1~ 2 |w| 3) oo \u\ 3 , 

and 

V 2 = t|u|i i0o + r 2 |w| 2iOC + T 3 \u\ 3yOQ + t 3/2 ||u|| Xt . 
Let r > 9/2 be fixed. The Sobolev embedding theorem, the two lemmas above, and ([7]) imply 

|lWt)lk T(t) <f(t)|| W (t)||r T(t) +C|| W (t)||^(l + r(t) 2 ) 
+ C\\u(t)\\ YT(t) (r(*)||Vu(t)|U- + (r(t) 2 + r(tf)\\u(t)\\ Hr + r(t) 3 / 2 \\u(t)\\ x ^ t) ) . (8) 

If r(i) decreases fast enough so that for all < t < T* we have 

f (*) + Cr(*)||V«(t)|Uo. + C{r{tf + r(tn\u(t)\\ Hr + Cr(t) 3 / 2 \\u(t)\\ Xr(t} < 0, (9) 
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then ((SJ implies that 

dt' 



d ;\\u(t)\\x T(t) <C\\u(t)\\U^ + r(On 



and therefore 

2 



||«(t)||x rW < ||«o||x tCO j + C T (o) / \\u{s)Y H rds = M(t), 

Jo 

for all < t < T», where C T (o) = 1 + t(0) 2 . Since r must be chosen to be a decreasing function, a sufficient 
condition for (jjjJJl to hold is that 

f(i)+CT(i)||V M (t)|| L ^+Cr(i) 3 / 2 (c; (0) || W (t)||^+Af(i)) <0, (10) 

where C' r , Q -, = t(O) 1 / 2 + r(0) 3 / 2 . For simplicity of the exposition we denote 

G{t) = cxp (c J \\Vu(s)\\ L oods\ , 

where the constant C > is taken sufficiently large so that ||u(t)|||fr < ||uo||#rG(£). It then follows that 
(jTUf is satisfied if we let 

r(t) = G(t)- 1 / 2 {^r 1/2 + C [ {C T{0) \\u{s)\\ H r +M(sj) Gisy'ds) . 
The lower bound (UJ on the radius of analyticity stated in Theorem 12.11 is then obtained by noting that 
r(0)- 1/2 +cj (g; (0) |Ks)||^ + M(s)) G^ds 

< r(O)- 1 / 2 +cj (c; (0) ||«o||*r + ||«o||jc t( o) + sC T( o)|| Uo ||^) da 

<G (l + i) 2 , (11) 

and therefore 

Go 



r(i)>G(t)- 1/2 - 



1 + t 

The last inequality in (fTTj) above gives the explicit dependence of Go on u . This concludes the a priori 
estimates that are used to prove Theorem l2.ll The proof can be made formal by considering an approximating 
solution vS n ' , n € N, proving the above estimates for vS n \ and then taking the limit as n — > oo. We omit 
these details. 

4. The commutator estimate 

Before we prove Lemma |3. II we state and prove two useful lemmas about multi-indexes, that will be used 
throughout in Sections 2] and [3] below. 



Lemma 4.1. We have 

/or a// a, /3 G Nq with (3 < a. 
Proof. Using (0J we have that 



a'\\ (a 3 



M a Mp X M-i p < ( |"! ) (12) 



and hence the left side of (|T2j) is bounded by 
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The lemma then follows from 

n\ I m\ I n + m 



for any n, m > such that n > i and m > j, which in turn we obtain by computing the coefficient in front 
of x x+ 3 in the binomial expansions of (1 + x) n (l + x) m and (1 + x) m+n . D 

The second lemma allows us to re-write certain double sums involving multi-indices. 

Lemma 4.2. Let {xaJagn 3 an d {z/a}a6N 3 ^ e rea ^ numbers. Then we have 

5Z X! X py<*-I3 = I XI X P I ( H 2/-Y J - ( 13 ) 

|a|=m|)9|=j ) ^<a \l/8|=j / \ hi ="»-.? / 

The proof of the above lemma is omitted: it consists of re-labeling of the terms on the left side of (|13l) . 
Now we proceed by proving the commutator estimate. 



Proof of Lemma \3.1\ We have 



where we denoted 



oo m 



m=3 j=l 



C ^=Un—W^ E M Q " ||^. V9«-^|U 2 . (14) 

1 ; ' |<*|=ro |/8|=j,/9<a V 7 

We now split the right side of the above equality into seven terms according to the values of m and j, and 
prove the following estimates. For low j, we claim 



y^ C m s < C f |«|i )0O |tt|3 + C'r|u|i ) oo||'u||y T , (15) 

m— 3 

oo 

Y, c m,2 < C|«| 2l oo|«|2 + Ct\u\ 2 ^\u\ 3 + Cr 2 |u| 2j00 ||-u|| n , (16) 



m— 3 

for intermediate j, we have 

oo [m/2] 



J2 T,CmJ<CT 3 ' 2 \\u\\ XT \\u\W r , (17) 

m= 6 j— 3 
oo m— 3 

£ E ^^^NIxJMk. (18) 

m=7j=[m/2] + l 



and for high j, 



J2 C ^ra-2 < CT 3 |u| 3 ,oo||w||y T , (19) 

771—5 
OO 

} J C m ,m-1 < C , r|u| 2 ,oo|w|3 + CV 2 |m| 2: oc|M|y t , (20) 

771—4 

OO 

X C-™,,m < C|u|i j00 |m| 3 -I- Cr|u|i !00 ||u||y T . (21) 

771—3 

Due to symmetry we shall only prove (fT5]) - (fT7l) and indicate the necessary modifications for (fl8 ]H [2l |) . 
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Proof of (|T5|) : The Holder inequality, (HIT) , and Lemma 14.11 imply that 

oo 

E C »M = E E (MpWdPuU-) (M a _p\\d a -PVu\\ L 2) M a M^M-l fj 

m=3 l«l=3 \j3\ = l,P<a 

3)! 



EE E (^II^IM (w^iifl-^vt*^^^ 

m=4 |a|=m |/3| = l,/3<a ^ ^ 



^ C E E (Mp\\d fj u\\ L ~) (M a ^ p \\d a -^u\\ L 2 

\ot\=3 \P\=l,P<a 

oo 

+ CV E E E O^ii^uM 

m=4 | a |=m |/3| = l,/3<a 

-rn— 4 N 

x M„_* 



The first sum on the far right side of (|22[) can be estimated by 

C|m|i,oo|Vm| 2 < C|w|i ]00 |u| 3 . 
Since m > 4, Lemma T4.2I implies that the second term on the far right side of ([22]) is bounded by 



°2^ ( m __ 3-)_m-4 

Cr|u| l!O0 V |Vu| m _ii- '—^ < Cr|tt|i,oo||«||*v, 

*-— [m — 3)! s 



m— 4 



concluding the proof of (|T5|) . 

Proof of (fT6|) : As in the proof of (fT5J) above, we have 



E C ™,2<C Y, E T m - 3 (A/ /3 ||^ u || L oo)(M a ^||a Q ^V U || L2 

m=3 |a|=3,4|,3|=2,/3<a 

xM a M^M~^ 

+ c f:i: e (^ii^iu-) (M a _ g |ifl°-^ U || £a {m ~^iy~ 6 

m=5|a|=m|/3|=2,/3<Q ^ ^ ' 

xtf ° M ^(;) (m _ % _ 3) / 2 (23) 

Using Lemma 14.21 the first sum on the right of (|2"3")l can be estimated from above by 

C|u| 2 ,oo|Vli|i + CV|u| 2 ,oo|Vu|2 < C|u| 2 ,oo|w|2 + Ct\ u\ 2 ,oo \u\ 3 . 

On the other hand, since s > 1, \/3\ =2, and \a\ — m > 5, we have by Lemma [O] that 



M a M p M a _ fj \^j j- _ 4)(m - — < ^ 2 j (m - 4)(m - 3) < C. 
By Lemma [OJ the second sum on the right of ([23]) is thus bounded by 

00 ( m _ A\ T m-5 

Cr 2 ]T \ukooNu\ m ^ [ > < Ct 2 \u\%oo\\u\\y t . 

This proves the desired estimate. 

Proof of (fTTf : We first observe that the Holder inequality and the Sobolev inequality give 

W&u • Vd a -' 3 u\\ L , < C||^||^ 4 ||A^ U ||^ 4 ||V9 a -^|U 2 . 
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Therefore we can bound the right hand side of (IT71) as follows 

,11,, 

(i-3) 



oo [m/2] oo [m/2] , _ 3 , 1/4 

m=6 j=3 m=6 j=3 \a\=m \0\=j,0<a V 



x 



M p \\d p Au\\ L 2 .. — Ma^H^-^VttlU 



where 



A*,/3, s = M a M^M- Y _ 



(./-1)!V V P " (ro-j-2)!' 

(j-3)p/ 4 (j-l)! 3s / 4 (m-j-2)! s 



o 



* Q -^W (m-3)! ;5 (m-j-2) 

By Lemma T4.ll we have that for to > 6 and 3 < j < [m/2] 



}J\j-lJ {m-j- 2)0" - l) s/4 (j - 2) s / 4 



<c|- 3 ^ s+1 ' 



j-iy ji+s/2- 

Since s > 1 the above chain of inequalities gives that A a ,p, s < C. Together with Lemma |4~2"1 and the discrete 
Holder inequality this shows that 

oo [m/2] oo [m/2] , j _ 3 . 1/4 , ] _ l s 3/4 

E E <w < ^ 3/2 E E (NiTTivO " Au ^ tttt) ' 

m=6 J=3 m=6 j=3 V U '' ' V U ' 

/, , (m-j- 2)t" 1 -3- 3 
The discrete Young and Holder inequalities then give 

oo [m/2] 

£ £c mj <CT 3 / 2 H|x T N| n , 

m— 6 j— 3 

concluding the proof of flTTf . 

To prove (HHJ) - (EH) we proceed as in the proofs of (rr5j) - (fTT)) above, with the roles of j and to — j reversed. 
Instead of estimating \\d°u ■ V<9 Q_/3 w||£2 with ||9' 3 u||.E,°o||d Q_ ' 3 Vu||£2 we instead bound 

||0VV0 a -*u|| La < ||0 /9 u|M|0 o "- /9 Vv||£~. 
We omit further details. This concludes the proof of Lemma |3~T1 □ 

5. The pressure estimate 

In the proof of the Lemma 13.21 we need to use the following higher regularity estimate on the solution of 
the Neumann problem associated to the Poisson equation for the half-space. 

Lemma 5.1. Assume that p is a smooth solution of the Neumann problem 

—Ap = v in £1, 

— — = on <9f2, 
on 

with v € C°°(J7). Then there is a universal constant C > such that 

\\d 3 d a p\\ L *<C J2 (***W«IU'> ( 24 ) 

/3'-a' = (2s,2i) 
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for any m > 1 and any multiindex a G Nq with \a\ = m and a 3 7^ 0. Additionally, if ct 3 > 2 then 

Wdidyu* < c Yl ( s ^)ll^iu 2 , (25) 

P'~a' = (2s+l,2t) 

\\d 2 d a P \\ L i < c J2 (***W«IU*> (26) 

/3'-a' = (2s,2t+l) 

where C > is a universal constant. 

We emphasize that the constant C in the above lemma is independent of a and m. In (|24[) we have are 
summing over the set 

{(3 <E Nl : \/3\ = m - 1, 3s, t € N such that /3' - a' = (2s, 2t)} 

and similar conventions are used in (|2^|) . (|26p . and throughout this section. 



Proof. In order to avoid repetition, we only prove (|24|) and indicate the necessary changes for (|25|) and (|26[) . 
Let A' = <9n + c?22 be the tangential Laplacian. Using induction on k € No we obtain the identity 

df^p = (-A') k+1 p £ ^'(-A') fe -^, 

and upon applying d 3 to the above equation 

cf + 3 p = d 3 (-A') k+1 p- Y,d^ +1 (-A') k ^v. 

3=0 

Therefore given \a\ = m, with 03 = 2fc + 1 > 1, we have 

k 

d 3 d a P = <9 2fe+2 <9"';p = (-A'f+^V + X^ 1 )^^ 1 ^^ 11 +5 22 ) fc - J 5 Q 'v, (27) 

i=o 

and if 03 = 2fc + 2 > 2, we have 

fc 
d 3 d a P = di k+3 d a ' P = d 3 {-A') k+1 d a ' P + J2(-l) k - 3+1 dl 3+1 (d n + <9 22 ) fe ^cPV (28) 

3=0 

Since n = (0,0, —1), the function g = (—A') k d a p satisfies the Neumann problem 

-Ag = {-A') k d a 'v in SI, 

^- = on an. 

on 
Using the classical ff 2 -regularity argument for the Neumann problem we then have 

||A' ff || L2 <C||(-A') fe 9 Q '^|| L2 , 
and 

\\d 3 A'g\\ L 2<C\\d 3 (-A') k d a 'v\\ L 2, 
for a positive universal constant C. Combining the above estimates with (|27|) . (f28]h and the identity 



(dn+d^rw^Y, ( m )dl s dl m - 2s w, 



s=0 
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we obtain 

k 

\\d 3 d a p\\ L2 < Cj2\\dl J (dn + d 22 ) k -id a 'v\\ L 2 

3=0 

< CY / T / ( k ~ J )\\d 2 1 s+ai d 2 2 k - 2 '- 2s+a2 d 2 3 3 v\\ L2 (29) 

j=0 s=0 ^ ' 

if a 3 = 2k + 1 > 1, and 

k 

\\d 3 d a p\\ L , < CY, \\d 2 3 J+1 (dn + d 22 ) k ^d a ' v\\ L , 

3=0 

^ C E E ( k ~ j ) \\dl s+ai dl k - 2j - 2s+a2 df +x v\\ L , (30) 

j=0 s=0 ^ ' ' 

if a 3 = 2fc + 2 > 2. To simplify ([Ml) above, let t = k - j - s > and /3 = (2s + ai,2k - 2j - 2s + a 2 ,2j) = 
(2s + ax, 2t + a 2 , a 3 — 1 — 2s — 2t) £ Ng. Since |a| = m and a 3 = 2k + 1, we have \(3\ = m — 1, and by 
re-indexing the sums, (f2Tjf can be re-written as 

\\d 3 d« P \\ L z<c y, ( s +*Whi^ 

s,teN ,|/3|=m-l ^ ' ' 
P'-a' = (2s,2t) 

The above estimate also holds for a 3 = 2k + 2 with the substitution (3 = (2s + ai, 2fc — 2j — 2s + 02, 2j + 1), 
thereby simplifying the upper bound (1301) , and concluding the proof of (|24[) . 
To prove (123)) we proceed as above and obtain 

\\did a p\\ L2 = \\d^ +l d?di k+2 P \\ L * 

< C^Y, ( k ~ J )\\d* 1+2s+1 d 2 2+2k - 23 - 2s dl 3 v\\ L . (31) 

j=0 s=0 ^ ' 

if a 3 = 2fc + 2 > 2, and 

\\d x d a P \\ L * = \\d? +1 d?df +3 p\\ L * 

^ C EE ( k ~ J )\\d? 1+2s+l d2> +2k - 2j - 2s d 2 3 j+1 v\\ L * (32) 

if a 3 = 2fc + 3 > 3. In 1^ we let t = k - j - s > and j3 = (ai + 2s + 1, a 2 + 2i, 2j) = (a x + 2s + 1, a 2 + 
2£, a 3 -2-2s-2£),sincea 3 = 2fc + 2 and |a| = m. Similarly in (^ we let (3 = (ai +2s + l,a 2 + 2i, 2 j + 1 ) = 
(«i +2s + l, a2 + 2i, a 3 — 2 — 2s — 2£), since a 3 =2fc + 3and \a\ — m. The above substitutions and re-indexing 
prove (|2l)l) . Upon permuting the first and second coordinates, this also proves (|2l)l) . D 

Remark 4. We note that Lemma [57X1 does not give an estimate for ||9i9 a p||i2 and ||<92<9 Q p|J£2 if a 3 = 1. In 
this case we note that the function g = d a p satisfies the Neumann problem 

-Ag = d a 'v in Q, 

^- = on on. 

on 
The classical _ff 2 -regularity argument then gives 

\\dld a P \\ L 2 = \\dld 3 d a 'p\\ L 2 < C\\d a 'v\\ L 2, 

and 

\\d 2 d a P \\ L 2 = \\d 2 d 3 d a ' P \\ L 2 < c\\d a 'v\\ L 2, 

for a positive universal constant C > 0. 
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We note that Lemma IQ1 is different from the classical higher regularity estimates (cf. [THl I3U I3S] ) for the 
Neumann problem in the fact that the constant C in (|2"4")) - (l2l))) does not increase with m. The dependence 
on m is encoded in the sums with binomial weights on the right side of (|24p - (|26|) . 

The following lemma shows that only a factor of m is lost in the above higher regularity estimates if each 
||9"u||i2 term is paired with a proper binomial weight. This explains the definition of the homogeneous 
Sobolev norms | • | m in (j3|). 

Lemma 5.2. There exists a positive universal constant C such that 

[A/2] [ft/2] + _ & __ (, + t\ (fh + ft\ 

/or any to > 3 and any multi-index (3 = (fji, 02,171 — 1 — /3i — ^2) € Nq. Additionally, if Pi > 1 we have 

[(ft-l)/2] [ft/2] . . . A . . 

e e( a+ *: 2 2 ;_t )C)^( a a*)- ,34) 

while if 02 > 1 we /ia?;e 

[ft/2] [(ft-l)/2] . 

Y^ V^ (Pi+ P2-2s-2t- 1\ (s + t\ (Pi+P2\ ,„.. 

§ S 1 A -2. J I J SC H A )• (35) 

where C is a universal constant. 

We note that in particular the constant C is independent of to and /3. 

Proof. Due to symmetry we only give the proof of (1331) . Estimates (134)) and (|35p are proven mutatis-mutandi. 
First we recall that given a, 7 € Nq, with 7 < a, we have 

a\ (\a 



■JJ VI7I 
Using the above inequality we get 

[ft/2] [ft/2] + ^ _ 2s _ 2 . , + . , 0i + ^ ■ I 

s & v a-& A • A a 

< [ ^ 2] [ ^ 2] /A + A — A (A + AN - 1 < ^ [ g ] ^ + ^ - 1 

s=0 t=o ^ Pi ~ s / V Pi / s=0 t=0 

The lemma is then proven if we find a constant C such that 

[ft/2] [ft/2] _i 

E E T <CXA + A). 

s=0 i=0 ^ ' 

Without loss of generality we may assume that /3i, /?2 > 4. We split the above sum into 

[ft/2] [ft/2] . _1 [ft/2] _i . 1 [ft/2] 1 1 

ee cr ^ + " + E f /s+i 

s=0 t=0 v ' t=0 v 7 v 7 s=0 v 

[ft/2] [ft/2] _x 

+ E E C't 

s=2 t=2 V 



It is clear that 



Ti+T 2 + T 3 . (36) 

Ti+T 2 <C{Pi+p 2 ). (37) 
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We estimate T 3 by appealing to the Stirling estimate (cf. [3H p. 200]) 

e 7/8 Vn(~) <n!<e v / n(-) • 
This implies 



s\t\ <e9/s J st I 



( s + t)\- v s + t(i + s /ty (i + t/ s y 

Thus we obtain 



[ft/2] [/3 2 /2] _ 



Since s > 2, the Binomial Theorem implies 



and by (f38j) we have 

[/3i/2] [ft/2] /[ft/2] \ / oo -. \ 

t 3 <cj: £^<c E 1 )[Ew)* c *- 

s=2 t=2 \ s=2 I \t=2 / 

Since 0i + fi% < m — 1, the above inequality, (|36p . and ([37]) complete the proof of the lemma. □ 

Proof of Lemma \3. 2\ First, note that since p satisfies the elliptic Neumann problem (|P.1|) - (|P.2[) we may use 
Lemma 15.11 to estimate higher derivatives of d 3 p as 

J2 M a \\d 3 d a P \\ L 2<c J2 E M a f s+ s t ')\\d^d l u k d k u l )\\ L 2. 

\a\=m,a 3 ^0 \a\=m, a 3 ^0 s,t£N ,\P\=m-l ^ ' ' 

0'-a'=(2s,2t) 

By re-indexing the terms in the parenthesis, the right side of the above inequality may be re-written as 

[ft/2] [/3 2 /2] 

£ £ £ (* + J I " t ii^(a-*»Jiu.. 



101=01-1 s=0 4=0 v ^ 

Using the estimate |55|) of Lemma 15.21 we bound the above expression by 

Cm Y, MpWd^diUkdhudWv 

|/3|=m-l 

and therefore 



E E M a \\d 3 d" P \\ L2 * 



m— 3 



, (m-3)! s 

=3 \Ja|=m,a 3 ^0 ; 

^ C E E A^ll^(^fc^"«)ll^ r ^" 3 " ■ (39) 

m=3|/3|=m-l ^™ '' 

On the other hand, higher derivatives of d\p are estimated using the decomposition 

J2 M a \\did a p\\ L * = J2 M4did a p\\ L 2+ J2 MaWdtdPpWv. (40) 

|ck|— 772,0:37^0 |a|— rn,a3 — 1 |a|— 771,03 > 2 

By Remark [H the first term on the right of (HO)) is bounded by 

C Yl M a \\d a ' {d.UkdkuM^ = C Y, Mp\\d (dtuAu^W^. (41) 

M=m,ois=l \p\=m-l,p 3 =0 
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Using estimate (l2l>|) . the second term on the right side of (|4*u]) is estimated by 

/ \ 



|a|— m,a3>2 



i s,t&io,\p\=m—l 
\/3'-a' = (2s+l,2t) 



S+t 



J2 Ma )\\d^ (diU k d k Ui)\\ L , 



I 



By re-indexing the above expression equals 

[(/3-l)/2] [/3 2 /2] „,/,,, 

c e e E( A "i + _i;.r (r^ 1 ""'""" - 

|/3|=to-1,/3i>1 s=0 t=0 V A^l - / \ 

and using (IM)) it is bounded from above by 



Cm Y. MpWd? (d t u k d kUl )\\ L 2 

|y8|=m-l,^i>l 



Therefore, by (gOJ), (|41|). and (|42|). we have 

OO / 

J2 E Mjd^pWv 



(m-3)! s 



-3 \ |a|— m,a3^0 

oo 

<cj2 S M /9ii^( a *«fc^«<)iU a 

m=3 |/3|=m-l 



_m — 3 



(m-3)! s 



By symmetry, we also get 

£ I Y, M a \\d 2 d a P \\ L 2 



_m— 3 



(m-3)! s 



rn=3 \J Q |=m,a 3 ^0 

oo 

<C^ J] Mp\\d l3 {d l u k d k u l )\\ L 2 



r m — 3 



m=3 |/3|=m-l 

Combining (|39|) . (|43|). (|44|) . and the Leibniz rule we obtain 



(m-3)! s 



m=3 |0|=m-l 



,-m — 3 



(m - 3)!' 



oo m — 1 
m=3 j— 



where 



^-(Si 1 ^ ? ^(Qii^ft^-^^^ii 



L 2 - 
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(42) 



(43) 



(44) 



(45) 



|/3|=m-l| 7 |=j, 7 </3 

We split the right side of (|45p into seven terms according to the values of m and j . For low j , we claim 



^ ^m.o < C|w|i,oo|w|3 + CtIuIi^IIuIIy^ 

m— 3 

OO 

5^ Pm,l < C|w| 2 ,oo|u| 2 + Ct|w| 2 ,oo|u| 3 + Cr 2 |M| 2 ,oo||w||y T 

m— 3 

oo 

5^ ^m,2 < Cr 2 |u| 3 ,oo|w|3 + CT 3 |tt| 3 ,oo||«||y T 



(46) 
(47) 
(48) 



?n— 5 
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for intermediate j, we have 

oo [m/2]-l 
E E ■Pm, j <CT 3 / 2 \\u\\ XT \\u\\ YT (49) 

m^8 j"=3 

oo m — 3 

£ E ?Vi < C^IMIxJMk (50) 

m=6 j'=[m/2] 

and for high j, we claim 

oo 

E ^ni,m-2 < C , t|u| 2: oo|w|3 + CV 2 |u| 2 ,oo ||u|| Y T (51) 

m=4 
oo 

E ^m,m-l < CMi.ooMa + CT|u|l,ool|t*||lV- (52) 

m— 3 

The above estimates are proven similarly to (|15l) - (|21[) in the proof of Lemma |3. II Due to symmetry we have 
presented there the proofs of the estimates where j < m — j. For completeness of the exposition we provide 
the proofs of (|5(T|) (|5^|) , where we have m — j < j. 

Proof of ([50)1 : We proceed as in the proof of (fTT|) in Section 0J First, the Holder and Sobolev inequalities 
imply that 

\\andiUk ■ df-idkUiWv < cwdW^w^Wd^dku^WAd^dku^ 4 . 

Therefore, 

oo 771—3 oo m— 3 / f ' n\ j — 3\ 

E E ^^E E E E [M^oM\r, { ^%A 



m=6 j=[m/2] m=6 j=[m/2] |/3|=m-l |7|=i,7</3 



x Nkf^ll^-T^H 



r m-j-3 \ 1/4 



L 2 



(m-j-3)!« 



where 



—rn—j — l \ 3/4 



R M M-iM-i ^ m(j - 2)! s (m - j - 3)!*/ 4 (m - j - l)! 3s / 4 



By Lemma |4. II we have that for m > 6 and [to/ 2] < j < to — 3 



S/3, 7 , s < C 



— s 



j J\j-2j (j - 2)(to - j - l) s / 4 (m - j - 2) s / 4 

<c(™r 2 3 ) 1_ V-:/r* /2 , 

since ( m ~ ) < C^J^o)) when j > to/2. Therefore, Bp,-y, s < Cj hence, by Lemma 1431 and the discrete Holder 
inequality, we have 

oo m — 3 oo m — 3 / m— 7 — 3 \ 1/4 

£ £ V m ,<Cr^Y. E (l^U-J-i (m T _.-_ 3) , 

ro=6j=[m/2] m=6j=[m/2] V V ; ' 



r 



m-i-l \ 3 / 4 



X |A9fc«i|m-3-17 : 7T77 |<9iW/c 



(m-j-l)!^ V (i- 2 ) ! 



(.7 - 2)r^ 3 
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The discrete Young and Holder inequalities then give 

oo m— 3 

£ £ ■p m , J <CT 3 / 2 \\u\\ Xr \\u\\ Yr , 

m=6 j=[m/2] 

concluding the proof of (fSOT) . 

Proof of (15T1) : As above we use the Holder inequality and obtain 



f>™,™-2<E E E MJ^wo^diUkW^wd^dkUi 

ro=4 m=4| / 9|=m-l|7|=m-2,7</3 ^'' 

^ Ct E E ll9 7 a^|| L2 ||9^9 fcUl || L ^ 

|,3|=3 | 7 |=2, 7 </3 

+ ^ 2 E E E (^ii^^-ii^T^i 

m=5 |/3|=m-l |7|=m-2,7<y3 V *■ ; 

x (M^ 7 \\d^d kUi \\ L oo) MpM- x M^ (P > ' 



(to-3)! s 



^7/ (to — 3) s 

Using Lemma FTTl Lemma [OJ and s > 1, this shows that the far right side of the above chain of inequalities 
is bounded by 

00 ™_m— 5 

Cr|9iUfc| 2 |9feUi|i oo + Cr 2 |5fcM l |i, 0O Y^ |9jU fe | m _ 2 7 7777 

*— ' (to — 4)! s 

m=5 v ; 

< C , t|u| 2 ,oo|u|3 + Ct 2 |u|2,oo||m||f t , 

thereby proving (f5T1) . 

Proof of ([5^)1 : By the Holder inequality we have 

00 oc m — 3 

m=3 m =3|y3| =m _l ^ m '" 

00 m — 4 

< C|9iUfc| 2 ||3A;Ui|| L oo + Cr ||<9fcUi|| L oo 2J |^iWfc|m-i7 — _ „x |s 

m=4 "• '" 

< C'|u|i i oo|u|3 + C , t|u|i i00 ||u|| Yt , 

which gives the desired estimate. By symmetry, we may similarly prove (|46[) - (|49p . but in these cases we 
apply the Holder inequality as 

\\d" f d i u k -d ~" f d k u l \\ L 2 < ||<9 7 a i u /£ || L ^||<9' 9 ~ 7 <9 fc u l || L 2, 

that is we reverse the roles of j and to— j. We omit further details. This concludes the proof of Lemma l3.2l □ 
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